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We prove Levinson's theorem for scattering on an (m + n)-vertex graph with n semi-infinite paths each 
attached to a different vertex, generalizing a previous result for the case n = 1. This theorem counts the 
number of bound states in terms of the winding of the determinant of the S-matrix. We also provide a proof 
that the bound states and incoming scattering states of the Hamiltonian together form a complete basis for 
the Hilbert space, generalizing another result for the case n = 1. 

I. INTRODUCTION 

<N 

t-H Continuous-time quantum walks on graphs were introduced by Farhi and Gutmann as a framework for developing 

new quantum algorithms. 6 This idea was subsequently applied to give an example of exponential speedup by quantum 
walk 3 and an optimal quantum algorithm for evaluating game trees. 5 Recently, it was shown that even a highly 
^ restricted model of continuous-time quantum walk is universal for quantum computation. 2 

A key feature of the quantum walks considered in Refs. 2, 5, and 6 is that the dynamics can be understood using 
scattering theory. For certain infinite graphs, one can construct an analog of standard quantum scattering theory, 
defining an S-matrix as well as eigenstates of the Hamiltonian corresponding to scattering of a wave packet at some 
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in 



X 



momentum. 2 ' 6 



Reference 4 proved a version of Levinson's theorem for continuous-time quantum walks on graphs obtained by 

i- 1 attaching a semi-infinite path to a single vertex of a finite graph. As in Levinson's original work, 10 the theorem 

*£h proven in Ref. 4 gives a relation between the phase of the reflection coefficient and the number of bound states 

supported by the Hamiltonian. In this paper we prove a more general version of Levinson's theorem for graphs. 

We begin by reviewing scattering theory and continuous-time quantum walks on graphs. At the end of this section 

"ti we summarize our new results and their relationship to previous work. 



Quantum walk and scattering theory 

(N 

Quantum walk is a quantum mechanical analog of classical random walk. In this paper we consider continuous- 
time quantum walk, 6 in which time evolution of the quantum walker occurs via the Schrodinger equation with a 



time-independent Hamiltonian. 



w~) A simple example is the quantum walk on an infinite path. The Hilbert space is spanned by basis vectors 

{\x):x£Z} 
Cn and the Hamiltonian is 



H= £ (\x)(x + l\ + \x + l)(x\). (1.1) 



Starting in a particular vertex state \y) and evolving with the Hamiltonian H, the quantum state after time t is 



-iHt 



\v)- 



To understand the dynamics, consider the analogy with the Schrodinger equation for a free particle in one dimension. 
There the Hamiltonian is 

H p2 

-"free — ~ — 

2m 

where the momentum operator p can be written in the position basis as —ih4-. The Hamiltonian (1.1) is simply the 
finite difference approximation to -fff rC o, up to an overall constant and a term proportional to the identity. As in the 
case of the free particle, the Hamiltonian is diagonalized in the basis of momentum eigenstates 



\k)= J2 e"***!*}, *e[-7r,4 




FIG. 1. Quantum walk on an infinite path with an obstacle represented by the directed graph G. The graph G includes the 
two vertices that are attached to the semi-infinite paths and may also include some internal vertices (not pictured). 

These states satisfy 

(k\p) = 2irS{k-p) 

where 6 is the Dirac delta function. Such a state \k) is an eigenstate of the Hamiltonian with energy E — 2cos(fc). 
The interpretation of these states is analogous to the interpretation of the momentum states for the free particle in 
one dimension. We can imagine preparing a wave packet, that is to say, a superposition of momentum states with 
momenta close to some value k. For k G (— n, 0) the wave packet moves to the left under the Hamiltonian evolution 
and for k G (0, n) it moves to the right. Indeed, the group velocity of such a wave packet is 

dE „ • MN 

~dfc m( } - 

Note that, unlike the case of a free particle, the dispersion relation E(k) = 2cos(fc) has the property that its phase 
velocity —E/k does not always have the same sign as its group velocity. However, it is the group velocity that 
determines the overall motion of a wave packet that is narrowly peaked in momentum space. 

More interesting solutions to the one-dimensional Schrodinger equation arise when a potential term V(x) is added 
to the free Hamiltonian H[ YCC . In the familiar case where V{x) goes to zero sufficiently fast as x — > ±oo, the eigenstates 
of the Hamiltonian 

P 2 

are of two types: bound states and scattering states. The bound states are normalizable states with amplitudes that 
asymptotically go to zero as x — > ±oo, whereas the scattering states asymptotically approach momentum states in 
this limit. 

One can define a quantum walk on any graph G. The Hilbert space has basis vectors labeled by the vertices of G. 
If G is an undirected graph, we take the Hamiltonian Hq to be the adjacency matrix of the graph. More generally, we 
can consider quantum walks on directed graphs where each edge i — > j is weighted by a complex number Wij, subject 
to the constraint that Wji — u>*- for all i,j. Then define the Hamiltonian to be the Hermitian operator 



H G :=J2 w io\JM (1-2) 



To consider scattering on graphs, we construct a graph G by attaching semi-infinite paths to some of the vertices of 
a finite graph G. The associated Hamiltonian Hq (defined through equation (1.2)) is equal to the adjacency matrix 
of the graph on the semi-infinite paths, and is equal to H G within the original graph G. The case of two semi-infinite 
paths, depicted in Figure 1, is closely analogous to the case of scattering off a (finite-range) one-dimensional potential. 
As in the one-dimensional Schrodinger equation, eigenstates can be either scattering states or bound states. The 
scattering states are eigenstates of this Hamiltonian, but they can also be viewed as describing the dynamics of a 
wave packet that is prepared on one semi-infinite path and then allowed to evolve according to the Hamiltonian. Each 
scattering state has components that can be interpreted as an incident wave, a reflected wave, and a wave that is 
transmitted through the obstacle G. 

In this paper we discuss quantum walks on graphs obtained by attaching n semi-infinite paths to a graph G, as 
depicted in Figure 2. Now we can prepare an incoming wave packet on any of the n semi-infinite paths and allow it 
to scatter off of the obstacle G. Associated with each incoming momentum k G (— tv, 0) and each semi-infinite path 
j G {1,2,..., n}, there is a scattering eigenstate |sCj(fc)). The states {|sc,-(fc)) : j G {1,2,..., n}} can be compactly 
described by an n x n unitary matrix S called the S-matrix. 




FIG. 2. Graphs G that we consider in this paper consist of a finite, weighted graph G with m + n vertices attached to n 
semi-infinite paths. The m internal vertices of G are not pictured. We label the vertices on the semi-infinite paths (x,j) where 
» e {1,2,3,...} and je {1,2,. ..,n}. 

Outline of the paper 

In Section II we describe the bound and scattering eigenstates of the Hamiltonian, as well as the S-matrix, following 
Refs. 2 and 11. Theorem 1 of Section II (proved in the Appendix) shows that the incoming scattering states along 
with the bound states form a complete basis for the Hilbert space. As in previous work, 4,7 our proofs rely on analytic 
continuation of the S-matrix. We define this analytic continuation in Section III. In Section IV we prove Levinson's 
theorem for graphs of the form shown in Figure 2. Our proof relies on a technical lemma that is proven in Section V. 

Relation to previous work 

Levinson's theorem relates the number of bound states of the Hamiltonian to the winding number of det(S'(fc)) as 
k is varied from —n to n (i.e., the number of times the phase of the determinant wraps around the interval [0, 27r) as 
k varies). In the case n = 1, our theorem is simpler and slightly stronger than the theorem proven in Ref. 4, removing 
a minor technical requirement on the graph G. Our results can also be viewed as generalizing previous work which 
discusses discrete versions of Levinson's theorem on the half-line with a boundary condition at one end 1,8 (rather than 
a general weighted finite graph G as in Ref. 4). 

II. EIGENSTATES OF THE HAMILTONIAN 

Consider the Hamiltonian describing a quantum walk on a graph G as shown in Figure 2. Let G have m + n vertices, 
with n vertices attached to semi-infinite paths and m "internal" vertices. The Hamiltonian Hq associated with the 
graph G is defined as in equation (1.2). For convenience, we denote by H the (m + n) x (m + n) matrix of Hq in the 
basis of the m + n vertices of G. (In general, we use a hat to denote restriction to the finite graph) . We can write 

H ~\B D 

where the matrix A is n x n, B is m x n, and D is m x m. The full Hamiltonian H associated with the graph G is 
Hq plus a term connecting adjacent vertices on the semi-infinite paths: 

n oo 

H = h g +J2J2 (i^x* + ui + \ x + uxm) ■ 



Bound states 

Bound states are the normalizable eigenstates of the Hamiltonian, and therefore have amplitudes on each semi- 
infinite path that go to zero as the distance along the line increases. Each bound state \<fi) has the form 

(x,j\4>) = a j z x - 1 (11.1) 

on the semi-infinite paths (for j £ {1, . . . ,n} and x £ {1, 2, 3, . . .}) for some ctj £ C and z £ (— 1, 1)\{0} (the boundary 
cases where z = ±1 will be discussed later). Any state of this form automatically satisfies 

{x,j\H\4>)= (z+-\{x,j\4>) 

for x > 2 and j £ {1, . . . , n}. A bound state must also satisfy the eigenvalue equation at the m + n vertices of the 
graph. Let a be a column vector with entries etj for j = 1, . . . , n and let j3 be a column vector of the m amplitudes 
of the state \<j>) that are inside the graph and not on the semi-infinite paths. Then the eigenvalue equation can be 
written as 

S5)GM!)-H)C 

We write this in a more compact form using the operator 

7M:=(=V,n: B L,). (I") 



zB zD 



Then the eigenvalue equation is 



*>(;- 

Given any (m + n)-dimensional vector \v) = I g I satisfying 

j(z)\v)=0 

for some z £ (— 1, 1), there is an associated normalized bound state \4> v ) defined through 

(x,j\4>v) = N v a j z x ~ 1 for x £ {1,2, .. .} and j £ {1,2, ... ,n} 
(iu|</>„) = N v j3 w for internal vertices of the graph w G {1, . . . , rr, 

The normalizing constant N v is 

^112 



Confined and unconfined bound states 

Some bound states may have zero amplitude on each of the n semi-infinite paths, corresponding to a = in the 
previous discussion. These have been called confined bound states 4 or bound states of the second kind 11 . A confined 
bound state \ip c ) has nonzero amplitude only on the m internal vertices of the graph that are not on the semi-infinite 
paths. Writing 

'Ml 

for the restriction of \ip c ) to the vertices of the graph, we see that this vector must satisfy the eigenvalue equation 

'A B^ 



B D 



He) = A C |V> C ) 



for some A c G R. This implies D/3 C = A C/ 5 C and ifi /3 C = 0. Define the projector onto the semi-infinite paths 

oo n 
x—l q—1 

and its restriction to the graph 

'1 



p 

" '~ ' 

which projects onto the n vertices connected to the semi-infinite paths. A confined bound state \ip c ) satisfies P n \4>c} = 
and 

1 (z)$ c ) = (z\ c -z 2 -l)$ c ). (II.4) 

The roots of the polynomial z\ c — z 2 — 1 are z c and z~ x where 

Zc = \ (A c - VW^) ■ 

If I A c I > 2 then both roots are real, with one inside the unit circle and one outside. The vector space spanned by the 
confined bound states with energies |A C | > 2 is therefore 

C> := span Utp) : (i/)\P n \i>) = and j(z)\fy = for some z with \z\ < l| . (II.5) 

If | A c | < 2 then both roots are on the unit circle and are complex conjugates of one another. The vector space spanned 
by the confined bound states with energies |A C | < 2 is 

C< := span[|^): (ip\P n \i/)) = and 7(3) |$) = for some z with \z\ = 1 and z £ {-1,1}} . (II.6) 

If A c = ±2 then there is one repeated root equal to ±1. The corresponding vector space spanned by the confined 
bound states is 

C= :=span{|V>): (M p n\ip) =0 and 7 (z)|^) = for z e {-1,1}}. (II.7) 

In general, the Hamiltonian can have confined bound states as well as bound states with nonzero amplitudes on 
the semi-infinite paths. We can always form an orthonormal basis of the bound states of the Hamiltonian consisting 
of confined bound states in C := C< ®C = ®C > and bound states in its orthogonal complement C . We refer to bound 
states in C^ as "unconfined bound states." 

Half-bound states 

Half-bound states 8 are unnormalizable states that are "almost" bound states. They are eigenstates of the Hamil- 
tonian taking the form (II. 1) on the semi-infinite lines with z G { — 1, 1} (where a 7^ 0). 

Scattering states and the S-matrix 

For each k £ (— it, 0) we define a set of n incoming scattering states {\sCj(k)) : j £ {1, . . . , n}} which have the form 

(x,j'\Bc j (k))=e- ik *6 j > J + e ik *S J ,j(e ik ) 

— Zi^i'i + z x S 3 'j(z) 

on the semi-infinite lines, where z = e lk . Such a state has energy z + - = 2cosfe. (Here z is on the unit circle; later 
we discuss the analytic continuation of the S-matrix to other values of z.) The label j indicates the semi-infinite line 
on which the state is incoming and x — 1,2,3, .. . indexes the distance along this line (with \l,j) corresponding to 
the vertex where the jth line connects to the (m + n)-vertex graph) . 



We now write the eigenvalue equations that determine the amplitudes of |sCj(fc)) at internal vertices of the graph. 
For each j, write these m amplitudes as a column vector ipj(z), and collect these column vectors into a matrix 

*(z) := {Mz) Mz) ■■■ Mz)). 

The following matrix equation determines \fr(z) as well as the scattering matrix S(z): 

A Bt\ H z -i + zS (z))\ + /(z- 2 + z*S(z))\ = f z+ l \ (V* + *S(*)) 



The lower part of this equation says 



*(*) = 1 {z- l B + zBS(*)) , 

£ + z — D 



which determines each ipj( z ) m terms of the scattering matrix S(z). The upper part determines the scattering matrix. 
We find 

A (z- 1 + zS{z)) + B^(z) + (z- 2 + z 2 S(z)) = ( z +l\ (z- 1 + zS{z)) , 

which gives 

S{z) = -Q{z)- 1 Q{z- 1 ) (II.8) 

where 

Q(z):=l-z^ + Bt_I__^. 

Recall that (for now) we are restricting our attention to values of z on the unit circle. For such values, the S-matrix 
is unitary. To see this, note that Q(z^ 1 ) = Q{z)^ and [Q(z), Q{z~ x )\ = 0, so 

5(^t = -Q( 2 - 1 )t(Q( 2 )-i) t 
= -Q(z)Q(z- 1 )- 1 
= -Q{z- 1 y 1 Q(z) 
= S(z~ 1 ). (II.9) 

We can see from equation (II. 8) that 5(z _1 ) = iS'(z) -1 . This establishes unitarity of S(z) on the unit circle. 

A complete basis of eigenstates 

The bound states and incoming scattering states form a complete basis for the entire infinite-dimensional Hilbert 
space. Let {\tp c ) '■ c = 1, . . . ,n c } be an orthonormal basis of the confined bound state subspace C such that for each 
c, P n \i>c) = and 

B D ) ^ = Ac '^ c )- 

Furthermore, let {\4>b) '■ b = 1, . . . , rib} C C be a basis for the unconfined bound state subspace. In Appendix VI we 
prove the following theorem, generalizing the n = 1 case: 7 

Theorem 1. Let v and w be any two vertices of the graph G. Then 

ii »b He \ 

w-\soj(k)){soj(k)\ +J2 \4>b){<t>b\ +^lV'c)(V'c| \w) = 6 VW . 

" 7r 6=1 c=l / 



III. ANALYTIC CONTINUATION OF THE S-MATRIX 

We would like to analytically continue the scattering matrix S(z) from the unit circle to the rest of the complex 
plane. To do this, we rewrite equation (II. 8) so that each matrix element of S(e tk ) is manifestly a rational function 
of z = e lk . Recall the definition of j(z) in equation (II. 2). As can be verified by direct calculation, we have 

f v (1 z& \( -Q(z) 
^= zD-z 2 - 



zD-z 2 -l 



B 1 



and 

j(z)- 1 = 
From these expressions we obtain 



~Q{zy l 

jn^rBQ(z)- 1 1 




-^r 1 i(D = (§ll A). (in.!. 



'(*) -h 



Now consider this equation for z £ C (no longer restricting to the unit circle). Each matrix element of j(z) is a 
polynomial in z. It is invertible everywhere in C except at a set of points determined by the roots of the polynomial 
det(7(z)). We can write its inverse as 

det(7(z)) 

where adj(7(z)) is the adjugate matrix of "f(z). The matrix elements of adj(7(z)) are polynomials in z. Hence the 
entries of 7(z)~ 1 are rational functions of z, and so are the entries of — 7(z) -1 7(-). We therefore define the analytic 
continuation S(z) through equation (III. 1) , so S(z) is the upper left n x n submatrix of — 7(z) _1 7(-). 

Lemma 2. Let S(z) be defined through equation (III.l). Then each matrix element of S(z) is a rational function of 
z. Furthermore, if Zq is a pole of some matrix element of S(z), then either zq — or det(7(zo)) = 0. 

Proof. It follows from the above discussion that each matrix element of S(z) is a rational function of z. By equation 
(II. 2), the entries of the matrix z 2 j(-) are polynomials in z. Hence each matrix element of S(z) is a rational function 
of z with denominator z 2 det(7(z)). □ 

IV. LEVINSON'S THEOREM 

Levinson's theorem counts the number of bound states. Let us now make this more precise. Using expresions (II. 5), 
(II. 6), and (II.7), we define the number of confined bound states as 

n c := dimC = dimC> + dimC< + dimC = . 

Furthermore, define the number of unconfincd bound states as 



n b 



y^ dim (span j \ip) : \ip) € C> and ~f(x)\ij)) = o|j 



x£(-l,l) 

dct7(a;)=0 



and the number of (unconfincd) half-bound states as 



Y^ dim (span j | tf>) : \i/)) € C± and j(x)$) = o} 



n h := 

x=±l 

For our purposes, the half-bound states are only counted as half a bound state each. In other words, we consider the 
"number of bound states" to be 

1 
n c + n b + -n h . 



We now give another formula for the number of bound states. Define 

W{z) := dct( 7 (z)). 

Note that W(z) is a polynomial in z. Denote the multiset of roots of W by 

{zi,...,z k } 

where k is the degree of W (each root appears in the above list a number of times equal to its multiplicity) . Let 

ai ■= \{i: \zi\ < 1}\ 
2a 2 := \{i: \ Zi \ = 1 and z, £ {-1, 1}}| . 
a 3 :=|{*:zie{-l,l}}|. 

The following lemma relates the number of bound states to «i, a 2 , and a^. 
Lemma 3. With the definitions given above, 

Oil = Tib + dimC> 

a 2 = dimC< 

a 3 = n h + 2dimC = , 

so a\ + a 2 + |a 3 = nb + n c + \rih- 

The proof of Lemma 3 is given in Section V. 

We recall some useful facts from complex analysis. Given a closed, positively-oriented curve k in the complex plane 
and a complex function /(z) that is meromorphic in C and has no zeros or poles on k, we define the winding number 
w K (f) of / around k to be the number of times the image of k wraps around the origin. In other words, it is the 
number of times the complex phase of / wraps around the interval [0, 27r). The argument principle is a formula 
relating w K (f) to the number of zeros and poles of / inside the contour. It says that 

w K (f) = Z K (f) - P K (f) 

where Z K (f) is the number of zeros of / inside the contour n and P K (f) is the number of poles of / inside k (both 
counted with multiplicity). We also use the notation ^ K \{ a }(/) (respectively, P K \{ a }{f)) to indicate the number of 
zeros (respectively, poles) of / inside k but excluding the point a. 
We now prove Levinson's theorem: 

Theorem. The winding number of the determinant of the S-matrix around the unit circle T is 

wr(det(Sj) —2\m — n\ > — n c ~ -n/, 

This generalizes the main result of Ref. 4, which establishes the n = 1 case. Note that the approach of Ref. 4 has 
a technical requirement on G, namely that A ^ or B^B ^ 1 (in this case A is 1 x 1 and B is m x 1). Our theorem 
has no such technical requirement, regardless of the value of n. 

Also note that the form of Levinson's theorem depends on the conventions outlined in Section II for the definition 
of the S-matrix. A different convention for the relative phases of the scattering states would modify the statement of 
the theorem. 

Proof. Using equation (III. 1) , we find 

det(S(z)) = {-l) n z 



njlm " \z> 



W{z) 
= (-l) n z 2m 

= (-l) n z 2m 

= (-!)"( I ! :;| :'-'"'"'\;,. ' V " -■ <IV.1> 




9 

(Note that \ Zi \ > for all * because 7 (0) = -1, so W{0) = (-l) m+n ^ 0.) 

Although W may have roots on the unit circle, det(S(z)) does not have any zeros or poles on the unit circle. Indeed, 
|det(5)| = 1 on the unit circle since S is unitary there. We can also see this explicitly from equation (IV. 1). Note 
that each root Zj — ±1 is a root of both W(z) and W(-) and the 03 factors corresponding to these roots cancel in 
the ratio 

rij=i ( z - z j) 

appearing in the expression above. Similarly, for each root Zj such that \zj\ — 1 and Zj £ { — 1, 1}, there is another 
root z* = — . (For z £ R, j(z) is Hermitian and W(z) € M, so any roots of W with nonzero imaginary part must 

occur in complex conjugate pairs.) These 2a 2 roots also cancel in the expression for det(S*(z)). 

Similarly, if there are other roots Zj with < \zj\ < 1 such that — is also a root of W(z) then the corresponding 
factors in the numerator and denominator of equation (IV. 1) cancel each other. Let q be the number of such roots Zj 
(so that the total number of canceling factors in the ratio (IV. 2) is 2q). Then 

-Pruo} ( det S) =ai-q 

Z T \ {0} (det S) = (degree of numerator) — {jj= of zeros of numerator outside the unit circle) 

= (k-2q- 2a 2 - a 3 ) - («i - q) 

= k— q — ai — 2a 2 - a 3 , 

where "numerator" and "denominator" refer to the ratio in equation (IV. 2) after common factors have been canceled. 

The determinant of S may also have a zero or a pole at z = 0. From equation (IV. 1) we see that the lowest-order 
term in the Laurent expansion about z = is 2m — k, corresponding to a zero if this quantity is positive or a pole if 
it is negative. 

Finally, the argument principle shows that 

w T (det S) = Zr\{o} (det S) - Pr\{o} (det S) + 2m - k, 

— (k — q — a\ — 2a 2 — a 3 ) — (a% — q) + 2m — k 

= 2m — 2ot\ — 2«2 — 03 

1 
-n,b — n e - -n h 

where in the last line we have used Lemma 3. □ 

V. PROOF OF LEMMA 3 

The operator 

7W = (is l zD % _ x ) - AK -D + zB-i (v.i) 

has appeared in our discussions of the scattering and bound states of the Hamiltonian. In this section we first establish 
some technical properties of this operator that we use in the proof of Lemma 3. We use the following result of Kato 
(Theorem 6.1 and Section 6.2 of Ref. 9). 

Theorem 4. Suppose Tq,Ti,T 2 are N x N Hermitian matrices and consider 

T(x) = T + xT x + x 2 T 2 

as a function of the complex variable x. For each real x there exists an orthonormal set of eigenvectors {\wi(x)} : i = 
1,2, ... , N} ofT(x) which can be chosen to be holomorphic functions of x on the real axis. 

We now use this theorem to establish that the eigenvalues and eigenvectors of 7(2;) are smooth functions of x for 
x <G R. In fact, we show that one can choose a smooth basis for the eigenvectors that includes the confined bound 
states as basis vectors. In the proof of Lemma 3 we use this fact to write det(5'(x)) (for x € E) as a product of two 
terms: one term that incorporates the contribution of the confined bound states and another term that comes from 
the unconfined bound states. 



10 



Lemma 5. Let {\ipc}} be an orthonormal basis of confined bound states as described in Theorem 1. For x € K, 
j(x) is Hermitian and there is an orthonormal basis {\vi(x)) : i — 1, . . . , m + n — n c } and eigenvalues {ei(x) : i — 
1, . . . , m + n — n c } that are holomorphic functions of x on the real axis, such that 



m-\-n—n c 



70) =^Hc){i>c\{xXc-x 2 -l) + Y^ ei(x)\vi(x))(vi(x)\ forxeR. 

C=l 4=1 

Proof. Write 

7 («) = E i&x&i ( zX ° - z2 - J ) + h - E i^eX^oi 7(*o i - E I^X^cl ■ 

c=l V c=l / \ c=l / 

The second term can be written as 

M - JT $ e ) ($ c \ J 7 (z) f 1 - ]T $ c ) $,| j = M + zM 1 + z 2 M 2 

for Hermitian matrices Mo,Mi,M2 that are independent of z. The result follows by applying Theorem 4. 

We also use the following two lemmas in the proof of Lemma 3: 
Lemma 6. If e,(:Eo) = for some iq£R then 

1 



dej 
dx 



■'•() 



Xo + X (Vi(xo)\P n \vi(x )). 



Proof. Note that 7(0:0) \vi(xo)) = implies x n ^ (since 7(0) = —1). Dividing through by xq gives 

H\vi(x )) = — \vi(x )) + Xo ( 1 - P„ ) |^(xo)). 



We have 



so 



Now 



d/y 

dx 



d~j 
dx 



\vi(x )) 



H-2x(l-P n ), 



— -x (l-K) 
Xq \ J 



\vi(x Q )). 



dx e t (x) = {vi{x)\ I ^7(*)j !$(»)) + f ^(«K*)l] 7(a0|6i(a0) + (^'WbW f ^I^W) 



D 



(V.2) 



since \vi(x)) is a normalized eigenvector of 7(3;) (this is sometimes called the Hellmann-Feynman theorem). Using 
equation (V.2), 



dx 



Ci(x) 



1 

Xq 



Xo + Xo{Vi{xo)\P n \Vi{Xo)) ■ 



n 



Lemma 7. 

(a) IfW(z) = and < |z| < 1 z7ien Im(z) = 0. 

(b) IfW(z) = 0, |z| = 1, and Im(z) 7^ 0, t/ien (■0|P„|'0) = for any | - 0) to £^e null space 0/7(2;). 



11 



Proof. First consider part (a). It is clearly true when z = 0, so suppose \z\ > 0. The hypothesis W(z) = implies 
that there exists a normalized state \w) such that 



z 1 -y(z)\4>) = 



H - z 



so 



('/•I 



H 



1 



*P n 



^^n 



IV') = o, 



IV) = o. 



Writing z = re 1 ^ and taking the imaginary part of the above expression gives 



r(ip\P n W 



1 

r 



0. 



Since (^>|P„|'0) > 0, the bracketed expression is strictly positive when r € (0,1). Hence sin0 = 0, and therefore 
Im(z) = r sin ip = 0. 

For part (b), if \z\ = 1 but Im(z) ^ then the above equation (which holds for any \ip) in the null space of j(z)) 

says that ($^|$) =0. D 

We are now ready to give the proof of the main technical lemma. 

Proof of Lemma 3. We can use Lemma 5 to write 

/m-\-n—n c \ n c 



W(x) 



n ^) iik 



1) for x G 



(V.3) 



where ei(x) = (vi{x)\^(x)\vi(x)) are holomorphic functions of x for x G M. The above expression explicitly separates 
out the contribution of the confined bound states to the determinant. 

We now show that a\ = rif, + dimC>. Part (a) of Lemma 7 shows that all of the roots of W with magnitude less 
than 1 are real. Consider one such root Xq G (—1,1). If Xq is a root of the polynomial 



''A,. 



x 2 ~l 



then the other root is — ^ ( — 1, 1). So each confined bound state with energy A c = xq + is responsible for a zero 
of multiplicity one in the polynomial W(x) at x — xq, and no other zeros in the interval (—1, 1). Now turn to the 
unconfincd bound states. For each i — 1, . . . ,m + n — n c , using Lemma 6 we have 



dei 

dx 



1 
Xq 



- Xq ) + X (Vi(x )\Pn\vi{xo)) 



^0 



for xq G (—1,1) since the right-hand side has the same sign as xo (and xq ^ because 7(0) = —1). So each state 
\vi{x$)) such that e^Xo) = contributes a zero of multiplicity one at xq to the polynomial W(z). Therefore, for each 
xq G (—1, 1) such that W(xo) = 0, 

multiplicity of the zero at xq = dim I span < \ip) : \ip) G C > and 7(^0) W) = \ ) 



Hence 



ai 



E 


r e(-i,i) 


W(x )=0 



+ dim (span j|i/>): \ip) G C> and j(x )$) = Of) 
dim (span J I -0): \ip) € C> and j(x )\^) = o|) 



+ dim (spanj \tp) : \tp) G C> and j(x )\ip) = > 
rib + dim C> , 
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where in the last line we have used the fact that confined bound states corresponding to different energies are linearly 
independent (since D is Hermitian). 

We now show that a^ = n/> + 2dimC = . To understand the zeros of W(x) for x € { — 1, 1}, we again use equation 
(V.3). Each bound state in C = (i.e., each confined bound state with energy ±2) contributes a zero of order two located 
at x = ±1, since 

±2x - x 2 - 1 = - (x T I) 2 ■ 

On the other hand, there can also exist states \v~i(±l)) such that ei(±l) = 0. These are half-bound states and satisfy 
(by Lemma 6) 



dej 
civ 



= ±i<tr i (±i)|p n |tr i (±i)) 9 £o, 

±i 



so each such half-bound state contributes a zero of order 1. Hence 

a 3 = |{t:zi€{-l,l}}| 

= Y^ dim (span ||V^) : \ip) e C± and j(x)$) = o\) 
x=±\ <- 

+ 2dim (span l\i/>) : \ip) € C = and j(x)\^) = o\) 
= rih + 2 dim C— . 

Finally, we show that ot2 = dimC<. Now we are concerned with roots of W(z) on the unit circle but not on the 
real axis. We can no longer use the expression (V.3), so we now derive an alternate formula that can be used in this 
case. Let {|^ c )} be a basis of confined bound states as before. As in the proof of Lemma 5, write 

"c / n c \ / n c \ 

7 (z) = ]T|^)(^|(zA-z 2 -l)+ 1-£|&>$c| )l(z) 1-£|&><&| . 

c=l V c=\ ) \ c=l / 

From this expression we obtain 

n c 

W(z) = det(M(z)) IJ (zX c -z 2 -\) (V.4) 

c=l 

where M{z) is some (to + n — n c ) x (to + n — n c ) matrix. Let z$ satisfy \zq\ = 1 and zq ^ { — 1, 1}- From part (b) 
of Lemma 7 we see that all the corresponding bound states are confined bound states (i.e., satisfy (ip\P n \ip) = 0). 
This means that det(M(zo)) 7^ 0, so the number of zeros of W(z) at zq is the same as the number of zeros of the 
polynomial n"=i ( z ^ c — z 2 — l) at z . Each confined bound state with energy A c = z + i contributes a simple zero 
at zo and a simple zero at — = z'q 7^ zq, corresponding to the two roots of the polynomial 

zX c — z — 1. 

So the number of zeros of W(z) in {\z\ = 1 and z ^ { — 1, 1}} is twice the number of confined bound states in the 
subspace C<. In other words, 

2a 2 = \{i: \zi\ = 1 andz^{-l,l}}| 

= 2 Y dim ( s P an {W : WPn\ip) =0 and j(z)\4>) = oj) 

z: W(z)=0, 

\z\ = l, z^±l 

= 2dimC<. D 

VI. OPEN QUESTIONS 

Two directions for future work suggested previously 4 remain open. The first is to find an algorithmic use for 
Levinson's theorem in quantum computation. Another direction is to consider the inverse scattering problem, where 
the goal is to reconstruct as much information as possible about the obstacle G from scattered waves. 
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APPENDIX: COMPLETENESS OF SCATTERING AND BOUND STATES 

In this appendix we prove Theorem 1, establishing that the scattering states \sCj(k)) for j = 1, ...,n along with 
the unconfined and confined bound states form a complete basis for the Hilbert space. In other words, for any two 
vertices v and w in the graph, 



/0 i, «b «c \ 

s-MAOXsc^fc)! + E l&X&l + E^°W°I \w) = s m 

"* 6=1 o=l / 



This is a generalization of the n — - 1 case, 7 and our proof follows the same steps. The proof has three parts: 

1. Since vertices (r, q) and (s,w) with r, s > 2 (i.e., vertices outside the graph on the semi-infinite paths) have no 
overlap with confined bound states, i.e., 

(r,q\\T / \ip c )(iJ c \]\s,w) = 0, 
we prove 

(f° rlk H \ ( nb \ 

2. Similarly, for \r, q) with r > 2 and \v) with v a vertex in the (to + n)-vertex graph G, we show that 



/ Jj^ E Mfc))K(fc)| |«) = -(r,q\ f E 




3. When \v) and \w) are both basis states corresponding to vertices in G, the confined bound states play a role. 
We show that 



(w 



/^EK^KWi) \v) = S wv -(w\ ^El^^lj \v)-(w\ fEl^)^o|j |«). (A.2) 
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Part 1 



Here we consider vertices (r, q) and (s, w) outside the graph, so r, s £ {2, 3,4,.. .}. For q — w, we have 
r° dk 



r ^ /J^ w ' )Hb,w| |m) 



I ^ ( [e~ lkr + e lkr S qq (e lk )] [e*** + e _ifc *5 OT (e tt )*] + ^e ifc(r - s) 5 M -(e ifc )5 gj -(e ifc )* 



dfc 
.-,2?r 



(2 cos(A:(r - s)) + e ik ^ r+s ^ S qq {e ik ) + e~ lk ^+^S qq {e lk )*^ 



where we have used the fact that 5(e ) is unitary. Now, using the fact that 5(e) = 5(e )' (which can be seen 
from equation (II.9)), we get 



\r,q\ 



/•0 fjh. n \ r 

J ^I]l sc j( fc ))( sc j( fc )l \s,q)=5 rs + j 



71 dk 
27' 



,ik(r+s)g ( e ik ) 



(A.3) 



Now consider the matrix element (for q 7^ w) 

r° dk " 



., 5Zlscj(fc)><sCj-(fc)| |s,w) 

"■ 3=1 / 

/ 27 ( ^^ + e S i^ e )] e' ik8 S wq (e ik )* + e lkr S qw (e tk ) [e lks + e^ ks S ww (e lk ) 
J2 e ik ^- s ^S qj (e ik )S wj (e ik )A 



jg{q,w} 
dk 

^2n 



e - lk (r+s) Swg (e ik )* + e lk{r+s) S qw {e lk ) 



where again we have used the fact that S(e lk y = S(e tk ) l = 5(e lk ). Putting this together with equation (A.3), we 
have 



\r,q\ 



/ ^.^2\ SC j( k ))( 8C j( k )\ l S > W ) = S rs$qw + j 



v dk 

1 

2tt 



Mr+s) S qw {e lk ). 



(A.4) 



We now evaluate the second term on the right-hand side. Letting z = e , we write this as a contour integral over the 
unit circle T: 

<ik e* k ^S qw (e ik ) = I ^z r + s -'S qw (z)- 



2tt 



2iri 



Since S(z) is unitary on the unit circle, each of its columns and each of its rows is normalized to 1. So S qw (z) does 
not have any poles when z is on the unit circle. Furthermore, by Lemma 2, S qw (z) is a meromorphic function of z 
and all of its poles inside the unit circle occur at values z$ satisfying either z = or det(7(z )) = 0. By part (a) of 
Lemma 7 this means that all of its poles inside the unit circle lie on the real axis. Using the residue theorem, 

f ^i Zr+S ~ lSqw{z) = ^ (ReB ao [x r+ '- 1 S VB (x)]) 

residues xq&{— 1,1) 

where Kes Xo [f(x)\ denotes the residue of f(x) at x — xq. From (V.l), we have 
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using (III.l), this implies 



S pv {z) = (l,q\-y{z)- x ^)\l,w) 



(l,?l 



1 



-j+i-jKrft 



|l,u/> 



-4^, 



2; 



2 "<jto 



1- 



1 



(l,^^)- 1 !!,^}. 



(A.5) 



Thus 



dz 
2tti 



Z Jqw\Z) 



E 



Res, 



■** Uqw i~ •b 



r-\-s— 1 



l 



l 



(GKx)- 1 !!^) 



residues xoG( — 1,1) 

Since we are considering vertices outside the graph, r > 2 and s > 2. Thus the first term has no residues, so 



2?ri 



<-'gto( z ) 



E RCSa; o 



residues xq(E{ — 1,1) 



„r+s-l 



1- 



1 



{l.ffMar)- 1 !!,*;) 



We can use Lemma 5 to write 



m-\-n—n c 



2—1 V ' 



(A.6) 



c=l 



The confined bound states |-i/j c ) satisfy (ip c \l,j) = for all j E {1, . . . , n}. There is a residue for each state \vi(xo)) 
that satisfies "f(xo)\vi(xo)) — 0. We use Lemma 6 to evaluate these residues, giving 



dz 

2TTi 



Z £>qw\Z) 



E ReS ^o 



residues Xq€z(— 1,1) 



E E < 

zoG( — 1,1) i'- e;(a;o)=0 
VK(x o )=0 

E E - 

io£(- 1,1) «: ei(a:o)=0 
W(x )=0 



x r+s ~ x [ I 



7n+n— n c 



„r+s-l 



r+s — 1 



1 ^ ( l,g|^(o;o))(^(a: )|l,w) 

de^ I 
dx Ixo 



1 s 



(l,g|^(x ))(^(x )|l,w) 
x o/ ^ - x + xo{vi{x )\P n \vi{xo)) 



(A.7) 



As described in Section II, each vector \vi(xo)) such that ei(xa) = for some xq £ (—1,1) corresponds to an 
unconfined bound state \4>b) = \4>b(i,x )) an d vice versa. This related unconfined bound state agrees up to normalization 
with \vi(xo)) within the graph, so 

\<t>b(i,x )) = N i,x \vt(xo)); 
the amplitudes outside the graph are 

(y,j\<t>b(i,x )) = ^.^(i^l^iC^o))^ -1 for V e { 2 : 3 > 4 >---} and i e I 1 - 2 ,. ..,n}. 

From equation (II. 3), the normalizing constant is 



N h*o = i "2 + IK 1 - Pn)\Vi{x ))\\ 



l-xl 



1-xl 



1 -x% + xl(vi(x )\P n \vi(x )} / 
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where we have used the fact that (vi(xo)\vi(xo)) = 1. Substituting in equation (A. 7), we obtain 

dz 



jr-\-s — 1 



2iri 



Ko6(— 1,1) *: e;(zo)=0 
W(z o )=0 

= - E E a; O +S_2 ( 1 ! 9l^(^x o ))(06(i,a :o )|l,w) 

a 6(— 1,1) i: ei(x )=0 
W(x )=V 

= - E E ( r >#6(t,xo)X < /Vi,xo)h'" J ) 

xo£( — 1,1) t: ei(a;o)=0 
W(a:o)=0 

= -{r,q\\^\<f> b ){cj> b \\\ S ,w). 



Plugging this into equation (A. 4) gives equation (A.l), as claimed. 



Part 2 



We wish to evaluate 



dk 



* 27 Vi 



^(r,g|sc J (A:))(sc5-(fc)|z;) 



where \v) corresponds to one of the vertices in G, and where r > 2. 
We have 



71 

# - 2 cos(fc)) |s55 (fc)) = - Y, I CO (2, l\s Cj (k)) , 



i=i 



and 



<sc^(fc)|£>=E< sc ^ fe )l 2 >'><M 



-1 



1=1 



H -2cos(fc) 



J7„ « 



rfA- 



-1 



i=i 

/0 ,, n n 
^ E E O^* + e^^.(e^)) (e 2 *%, + e- 2 ^^(e^T) (1, l\ , 

-7T Z7T ~[ l ^ 1 H - 2 cos(/c) 

n 

v)+e^ r +^Y / S q i(e tk )(0\ 



v) 



* dk / <fc(r _ a)< Q 1 1 

. w 2tt I V ' y '#~2cos(fc) 



-1 



/=i 



H -2cos(fc) 



(A. 



where we have used unitarity of S as well as the fact that S(e lk ) j ' = S^e ,fc ). We can convert this expression to a 
contour integral over the unit circle, giving 

dk -A 

=r 

dz 



2tti \ V ,ql H-z- 



1 \v)+J2 z ^s ql {z){i;i\ 1 is)) 

z 1=1 z / 



2ni 



r-3 _ r-l 



M 



^r+1 



ff-*-i 



1 (1,317^)-^ 



F-z 



17 



where in the second line we have used equation (A. 5). Now by (V.l), 

Pn = \l{z)-- [H-Z-- 
Z A z \ z 



so 



Thus we have 



7(*)- 1 i > n - - 2 -l(z)- 1 (S z - ). 
z z z \ z) 



(A.9) 



J 77 'I' p J / -I \ 

w ^ £<r.9M*)><^-(*)|3) = / r S* r (* - ^J (QItW^R- (a.io) 

We now show that the integrand on the right-hand side has no poles on the unit circle, so that we can use the 
residue theorem to evaluate the contour integral. To see this, recall that 7(2) is block diagonal in an orthonormal 
basis that includes the confined bound states as basis vectors, where the two blocks correspond to the subspaces C 
(the subspace of all confined bound states) and C 1 - . Note that because (*0 C |1 5 <z) = 0, 

z r (1 - i) (Qmz)- 1 ® = 7/(1- 1) (Q| (1 - jr \j c )($ c \\ jiz)- 1 (1 - e |4><? C | J \v) 

Part (b) of Lemma 7 says that if 7(2:) has determinant at some z € T \ {±1} then its null space is a subspace of 
C. Hence the (m + n — n c ) x (m + n — n c ) block of j(z) restricted to C 1 - has nonzero determinant and is therefore 
invertible on T \ {±1}- Thus the integrand above has no poles on V \ {±1}. At the points z — ±1, it may be that 
det7(±l) = 0, but nevertheless 



1 

1 -z 2 



det ( --y(z) 



is nonzero for z £ {±1}- This follows from Lemma 6, which implies that an eigenvector in the null space of 7(±1) 
has a simple zero in its eigenvalue at z = ±1. So we have shown that the integrand in equation (A.10) has no poles 
on the unit circle. 

The poles of the integrand inside the unit circle occur on the real axis, where we can use the expression (A. 6) and 
Lemma 6 to evaluate their residues (as we did in part 1). This gives 

° dk^, I ,,«,-,Mn V- V- r/1 1 \ {Ql\vi{x )){vi(x )\v) 



[ SX>?M*)><^(*)|3> = E £ *o(l-i)x 

J ~' K 3=1 xoe(-l.l) »: e 4 (x )=0 V 0/ x 

W(x o )=0 



XO + x (vi (x ) J P n I Vi (xq ) ) 



= - E E X 1 ( 1 ><l\ ( f > b(i,xo))(<f>Hi,xo)\ v ) 

xoG(— 1,1) i: e;(x )=0 
W(x )=Q 

= -(r,«|(El^>(06|jl*) 



\b=l 



which completes the proof of part 2. 

Part 3 

To prove equation (A. 2) it is sufficient to show that 



"° dk n * / n " 



<x«l I / ^El sc ;( fc )>< sc i( fe )l \\x b ) = s ab -( X a\ (£l&X<fc,l 1 lx»>-{xa| (£hWc|J IX6> 
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for any orthonormal basis of m + n states {|Xa)} that are superpositions of the m + n basis states corresponding to 
vertices in the graph. We choose to work in the orthonormal basis of eigenstates of the matrix H, so 



A fit 
B D 



\Xa) =E a \x a ), 



and the state \x a ) in the extended Hilbert space is simply equal to \x a ) in the graph and has zero amplitude elsewhere. 
Therefore 

(r° dk n \ r° dk ™ 



and using equation (A. 8) we get 



<Ik f2(Xa\Sc j (k))(80 j (k)\X b )= f £ E 



2n 



i=i 



3,9,1=1 



(Xa\ 



1 



H -2cos(fc) 



1,«><M| 



1 



H -2cos(fc) 



IS) 



x ( e_2 *J« + e 2 S qj {e lk )) (e 2lk 5 ol + e~ 2lk : Sy(e lfe )*) 



dk 

2n 

4ik 



(Xa\ ft 777 P* ft 777 1 &) 

H-2cos(k) H-2cos{k) 



1 



e "^S <£a| ^-2cos( fc ) 



l,g)^(e lfc )(l,/| 



1 



27riz 



+ z 



(xol ^ r p »77 rl&>) 

H — z — - H — z — - 



1 



4 E<^I& T \l,q)S ql (z){l,l\ 



q,l = l 



H-z 



H-z 



H -2cos{k) 



\Xb) 



\Xb) 



By equation (III.l), 

Similarly to equation (A. 9), we have 

and 



S ql (z) = -(l,q\P nl (zrHl)Pn\lJ). 



p b ^)-i = ^-I(f-,-I'),u) l 



^^A(h-- rJ + TW'- 



r7O0 



3 - 1 I P - 



(A.11) 



1 1 

72 



F-^--) 7 ( z )" 1 



lilWn 



1 



z \ z 



:P n 



P n -UH-z-- z ) .(:. i 



r7(*) 



i^-ifi-iV^-^-ihw- 1 ^ 



-^-iri-iU^-,- 1 



1/1 



l\Pn 



p, 



-i)(H-z--) 1 (z)- 1 {H-z- 
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dk " 



where in the last step we have used equation (A. 9). Inserting this into equation (A. 11) gives 

n 

^(Xalscjik)) (aSj (k)\xb) 



2n 



dz 

r 2iri 



(A.12) 



Without loss of generality we assume the orthonormal basis {|Xa)} includes the confined bound states {|V> C }} since 
they are eigenstates of H . By equation (II. 4), 



i^r 1 ^) = 



i i 



z \ r — z - 



tI^c 



Plugging this into the above, we see that if either \x a ) or \\b) corresponds to a confined bound state then the integrand 
is zero. 

Now consider the case where both \x a ) and \x.b) are orthogonal to all confined bound states. Then the first term in 
equation (A.12) is 



dz 

r 2ni 



1 



1 



Sab 



E a -Z-\ 



* dk -2sm{k)8 ab 
_ 7T 2^E a -2cos{k) 



= 



since the integrand is an odd function of k. For the second term in equation 2, since the states \x a ) and \xb) are 
orthogonal to all confined bound states, the integrand 



z 1 



(Xah(z) 1 \Xb) = Z 



i - ^) (xa\ f i - f; m$c\\ i(z)- 1 f i - e i&x&i] 



x&; 



has no poles on the unit circle by the same argument given in part 2 for the integrand of equation (A. 10). As in 
parts 1 and 2, the second term in equation (A.12) has poles inside the unit circle where det7(z) = 0. However, 
unlike in parts 1 and 2, this term also has a pole at z — that must be taken into account. This pole has residue 
— (Xa|7(0) _1 |xb) = S a b, since 7(0) = — 1. Computing the residues at all of the other poles proceeds as in parts 1 and 
2, and we obtain 



£ £ K 1 _ ^) (£a|7(zri|£6> ] = Sab {xai (£ |06)(< ^) 



\Xb) 



Hence 



j;„ " 



dk 

2~i 



X](Xa|sC j (fc))(sC i (fc)| 



Xb) 



2=1 



if either |^ a ) or \xb) is a confined bound state 

Sab - (Xa\ (E?=i l&X&l) IXfa) if (XalVO = (ft|&) = for all c € {1, • • • ,n c }. 



Equivalently, 

(Xo| 



dk 

2^ 



E|sc J -(fc))(sc i (fc)| 



3=1 



\Xb) = S ab ( Xa \ £ |&><&| |X6) " (Xal £ IVcX^cl |X&>, 



\b=l 



which is what we wanted to prove. 



